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Abstract 

The AMS-02 collaboration has recently reported an excess of the cosmic- 
ray positron fraction, which turned out to be consistent with previous results 
reported by the PAMELA and Fermi-LAT collaborations. A decaying dark 
matter with the mass around 1 TeV can be responsible for the excess of the 
positron fraction when it is interpreted as a dark matter signal. Interestingly, 
the pure gravity mediation model provides such a dark matter, namely an 
almost pure neutral wino dark matter, when a tiny R-parity violation through 
LLE^ interactions is introduced. We show that the decaying wino dark matter 
well reproduces the energy spectrum of the fraction with being consistent with 
constraints from cosmic-ray anti-proton and gamma-ray observations. 



1 Introduction 



The AMS-02 collaboration recently released their first result of the cosmic-ray positron 
fraction [T]. They reconfirmed the anomalous excess of the positron fraction over 
the expectation based on simple cosmic-ray propagation models, which had been 
found by the PAMELA [2j and Fermi-LAT [3] collaborations. One interesting fact 
the AMS-02 collaboration found is that the positron fraction shows no observable 
anisotropy. Though it is obviously premature to make a definite statement on it 
because of limited statistics, this fact encourages several dark matter interpretations 
of the excess rather than astrophysical interpretations. 

When the excess of the positron fraction is interpreted as a dark matter signal, 
a decaying dark matter with the mass around 1 TeV could be responsible for the 
excess. In particular, the decay of a wino-like neutralino dark matter is known 
to naturally account for such an excess as discussed in reference j^. The wino- 
like dark matter is a natural prediction when the gaugino masses are dominated 
by the anomaly mediated supersymmetry (SUSY) breaking contributions [S]. Such 
scenarios have been widely discussed in the models of pure gravity mediation [6l [7], 
the models with strong moduli stabilization [8], the spread supersymmetry [9], and 
the minimal split supersymmetry [10], all of which explain the observed mass of 
about 126 GeV for the Higgs boson [TT]. Thanks to the relatively large gravitino 
mass m3/2 = 0(100-1000) TeV in these models, almost all of the phenomenological 
and cosmological problems in SUSY standard model are solved. The mass of the 
neutral wino dark matter is predicted to be m^ino = 0(0.1-1) TeV. In particular, 
the wino with the mass m„ino — 2.7 TeV is interesting since its thermal relic can 
explain the observed dark matter density if it is stable or long-lived [12] . 

In this paper, we show that the decaying wino dark matter with its mass around 
1 TeV can naturally account for the anomalous excess of the positron fraction re- 
ported by the AMS-02 collaboration. In particular, if the decay occurs dominantly 
through R-parity violating interactions LLE'^, the energy spectrum of the positron 
fraction is well reproduced with being consistent with constraints from cosmic-ray 
anti-proton and gamma-ray observations. After briefly reviewing the wino dark mat- 
ter in the framework of the pure gravity mediation model and discussing its decay 
processes through the LLE^ interactions in section [21 we explicitly show the param- 
eter region (the mass and the lifetime of the wino dark matter) which accounts for 
the AMS-02 result in section [31 Section [H is devoted to summary and discussions of 
some prospects for testing the decaying wino dark matter in (near) future. In the 
appendix A, we also discuss a "GUT" model for the R-parity violation. 
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2 Wino dark matter in the PGM model 



We first summarize the wino dark matter predicted by the pure gravity mediation 
(PGM) model [B] and its decay into electrons/positrons through LLE'^ interactions, 



#1 



In 



where L and E'^ are the lepton doublet and charged lepton singlet, respectively 
the PGM model, all gauginos acquire their masses from one-loop anomaly mediated 
contributions p], and the neutral wino becomes the LSP. The Higgsino mass term, 
i.e. /i-term, of order of the gravitino mass is generated through interactions to the 



R-symmetry breaking sector [13] |^ Thus, with such a large yU.-term, the mixing 
between the wino and the bino is highly suppressed. The PGM model therefore 
predicts the almost pure neutral wino as a dark matter candidate. 

The lower limit on the wino mass is currently set by the Large Hadron Collider 
experiment [15] and the Fermi-LAT [16] observation. In the former case, the direct 
charged wino production with a disappearing charged track [6l [17] gives a limit as 
'^wino ^110 GeV [TBj. In the latter case, the observation of gamma-rays from wino 
dark matter annihilations at milky-way satellites (dwarf spheroidals) gives a limit 
as 2.2 TeV > m„ino > 500 GeV and m„ino > 2.4 TeV [191. On the other hand, the 



upper limit on the wino mass is obtained by cosmology. The wino dark matter is 



produced thermally [12] and non-thermally through gravitino decays'^ i in the early 
universe, and both contribute to the dark matter density observed today. It then 
turns out that the wino mass must be lighter than ~2.7 TeV in order not to be 
over-produced. The wino mass of around 2.7 TeV is particularly interesting because 
the dark matter density is explained only by its thermal relic [T2] . 

The wino dark matter decays into standard model particles when the R-parity 
is violated in the PGM model. There are several ways to introduce the R-parity 
violation, which makes the wino meta-stable. The simplest way is the use of LiHu 
interactions with being the up-type Higgs doublet. Then the wino decays into 
a Z boson plus a neutrino, a Higgs boson plus a neutrino, and a W boson plus a 
charged lepton [21]. Though the decay processes indeed produce high-energy elec- 
trons/positrons, they also produce anti-protons because of cascade decays of Z, W 
and Higgs bosons, leading to a severe limit on its lifetime [22]. Furthermore, those 
cascade decays produce gamma-rays with a hard spectrum. High-energy charged lep- 



'^^The analysis developed in the paper is also applicable to the most of models with heavy sfermion 

as long as the lightest supersymmetric particle (LSP) is the almost pure wino-like neutralino. 
'^^This effect is not exactly the same as the original Giudice-Masicro mechanism [T4|, because the 

Higgsino mass term (the /Lt-term) is not generated by interactions to the SUSY breaking sector. 
#3The non-thermal contribution depends also on the reheating temperature (Tr) after inflation, 

where it is estimated to be n'^^^lh'^ ~ 0.16 x (m^ino/SOO GeV) x {Tr/W^ GeV) [20l|6]. 



2 




Figure 1: The decay of the wino dark matter through the interaction Xij^LiLjE'^. 

tons produced by the two-body decay processes also generate energetic gamma-rays 
through the inverse Compton scattering (ICS) with the cosmic microwave back- 
grounds and infra-photons from star hghts [23]. Those facts make the decaying wino 
dark matter interpretation of the positron excess difficult [23]. 

In this letter, we therefore focus on the wino decay caused by the interactions, 

>V^ = A,,,LiL,E^, (1) 

where we do not have to worry about the constraint from cosmic-ray anti-proton 
observations, while another constraint from gamma-ray observations also becomes 



milder than the case of LH^ [3]- 1^ The assumption that R-parity violation is 
dominated by LLE'^ interactions needs to be justified in the grand unified theory 
(GUT) in which quarks and leptons are unified into GUT multiplets. We discuss a 
model addressing the origins of the interactions in appendix |Al where only LLE'^ 
are generated while suppressing hadronic R-parity violating interactions in a GUT 
consistent way. 

Through the R-parity violating interactions in Eq. ([T]), the wino dark matter 
decays into two charged leptons by emitting a neutrino through three-body processes 



as shown in figure [H The wino lifetime is then estimated to be 

rwino ~ 10'' [sec] (A/10-i^)-2 {m^ino/l TeV)"^ (m^/lO' TeV)^ (2) 

The decay can be a source of high-energy cosmic-rays which results in the excess of 
the positron fraction. In the next section, we discuss positron, electron and gamma- 
ray fiuxes which are expected from the decaying wino dark matter. 



There also be some studies discussing contributions to anti-proton and gamma-ray fluxes from 
next leading order processes such as the electroweak bremsstrahlung, where these contributions 
become important in some models aiming to explain the AMS-02 result |25j . In our model, with 
the wino mass of 0(1) TeV, these processes are less significant [26] , 
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3 Signals of the decaying wino dark matter 



For dark matter signals, we assume that the left-handed slepton of the third gen- 
eration (I/3) is lighter than those of other generations (I/i,Z/2) and the wino dark 
matter decay occurs dominantly by exchanging L3. We also assume no flavor viola- 
tion on the coupling between Lj, Lj and for simplicity. In following discussions, 
we consider the cases of A32j [i = 1,2,3) as a demonstration. Branching fraction of 
the decay is then given by BF(iy° Tiu^ejn) = BF(iy° -)■ fiL^r^Ri) = 0.5. Sub- 
sequent decays of heavier leptons lead to energetic positrons, electrons, neutrinos 
and photons. High-energy photons are produced also from final state radiations. 
We estimate the primary spectra of the particles from these prompt decay channels 
using the Pythia 8 code [2B] with a slight modification by us to treat the decay of a 
polarized lepton. 

3.1 Electron and positron fluxes 

In the calculation of electron and positron fluxes from the wino decay, the NFW 
profile [21] is used for the dark matter mass density of our galaxy with the profile 
parameters being Pq = 0.4 GeV/cm^ (the local halo density), Vc = 20 kpc (the core 
radius), and Tq = 8.5 kpc (the distance between our solar system and the galactic 
center). In order to take account of the effect of electron and positron propagations 
inside our galaxy, we have solved the diffusion equations according to the method in 
reference [30]. We have mainly used the so-called MED model in the reference, but 
other models Ml and M2 have also been used to estimate the uncertainty associated 
with the propagations. The Ml (M2) model maximizes (minimizes) the fluxes with 
the Boron to Carbon ratio being consistent with its observation. 

For astrophysical backgrounds against the signals, we must be careful about their 
calculations because of large uncertainties. In order to deal with the uncertainties, 
we have adopted a similar method developed in reference [31]. Using four parameters 

and p^, the interstellar background fluxes are parameterized as 

<^tGiE) = A^E^^<l>:t,{E), (3) 

where $^^f(£') are reference background fluxes, and we have adopted conventional 
ones obtained by the GALPROP code [32] with the electron injection index being 
—2.66, which gives the best fit to the observed electron flux at PAMELA [33]. In 
order to take account of the effect of the solar modulation, which is important for 
low energy electrons and positrons, we have used the force-field method [31] in both 
signal and background calculations. 



4 



We have fitted above electron and positron fluxes to several experimental data. 
The data we have considered are the electron flux at PAMELA [33], the combined flux 
of electron and positron at Fermi-LAT [35], and the positron fraction at AMS-02 [1]. 
In the fitting to Fermi-LAT data, we have included the systematic uncertainty asso- 
ciated with the absolute energy scale S'pcrmi- As a result, the electron and positron 
fluxes are fitted to the data by varying the following nine parameters: the wino 
lifetime (rwmo), the wino mass (m^mo), the background parameters {A"^, p"^), the 
force-field potentials for electrons and positrons (0^), and the energy scale ((SFermi)- 
We have varied the parameters A"^, p^, and (f)"^ within the ranges of A^ G [0, oo], 

G [—0.5,0.5], and 0^ G [0,1] GV, respectively, while adopted the error of the 
scale Spermi to bc (JS'Fermi = ±0.1. All the errors associated with S'pcrmi and experi- 
mental data are assumed to be Gaussian and we symmetrize the asymmetric error, 
if necessary. The fitting has been performed in the energy range of -E > 5 GeV for 
AMS-02 positron fraction, E > 10 GeV for PAMELA electron flux, and E > 20 GeV 
for Fermi-LAT electron and positron total flux. 

3.2 Gamma-ray fluxes 

We have considered several types of gamma-ray fluxes from the decaying wino dark 
matter in order to investigate how severely the mass and width of the wino dark 
matter are limited. The fluxes we have calculated are the following: The isotropic 
component of the diffuse gamma-ray flux, the diffused gamma-ray flux from all sky 
including the region of galactic center but not including the region of galactic plane, 
gamma-rays from galactic clusters. We have also considered limits on the mass from 
the wino annihilation, which are also obtained from several gamma-ray observations. 
In what follows, we explain how the gamma-ray fluxes have been calculated and how 
we have put limits on the mass and width of the wino. 

3.2.1 Isotropic component of the diffuse 7-ray 

In order to calculate the isotropic component of the diffuse gamma-ray flux from the 
wino decay, both galactic and extra-galactic contributions have to be considered. For 
the galactic contribution, we have calculated the gamma-ray flux in the direction 
of the anti-galactic center according to the dark matter profile mentioned above. 
We have not included the contribution from the ICS in the galactic halo, because 
its estimation has 0(1) uncertainties associated with the diffusion of electrons and 
positrons. For extra-galactic component, we have calculated the contribution from 
the wino decay at all past redshifts, where both prompt decay and ICS with the 
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cosmic microwave background are considered. The calculated flux has been compared 
with the experimental data officially reported by the Fermi-LAT collaboration [36] , 
and we have put a limit on the mass and width of the wino by demanding that the 
signal flux from the wino decay does not exceed any Fermi-LAT data points by more 
than two sigma significance. 

3.2.2 All sky survey of the diffuse 7-ray 

For the all sky survey, the calculation of the diffused gamma-ray flux from the 
wino decay is essentially the same as that of the isotropic component. Only the 
difference is that the flux is now calculated at all directions. We have considered 
both galactic and extra-galactic contributions in the same manner mentioned above. 
The calculated flux has been compared with another experimental data reported 
again by the Fermi-LAT collaboration [37]. We have considered almost of all region 
of the sky including the region of galactic center but excluding other regions along 
the galactic plane. The strategy to put a limit on the mass and width of the dark 
matter is also the same as in the case of isotropic diffused gamma-rays. 

3.2.3 7-ray from galactic clusters 

We have focused on the Fornax cluster as a target region to detect the gamma-ray 
signal, because this cluster is known to give the most severe limit on decaying dark 
matters among several clusters [38l [39] . In the calculation of the signal flux, we have 
included the effect concerning the extent of the dark matter proffie, which results 
in a weaker limit than that using a point-source approximation. We have assumed 
the NFW proffie for the dark matter density inside the cluster and estimated the 
proffie parameters ps and with use of Virial concentration parameter (cvir) [^1^ . 
Virial over-density parameter (Avir) [42], and the cluster mass [l3l HI] according to 
reference [12|. Though 0(1) uncertainty exists in the estimation, we have adopted the 
following values for the profile parameters: Ps = 0.09 GeV/cm^ and = 0.14 Mpc. 
In this setup, the so-called astrophysical J-factor for the signal flux with a radius of 
one degree centered at the target is J dQ = 1.9 x lO^^GeVcm"^. 

In order to estimate the background from the diffuse gamma-ray, we have used 
iso_p7v6source.txt for its isothermal component and gal_2yearp7v6_v0 . f its for 
its galactic component, where both are provided by the Fermi-LAT collaboration |45j . 
For the background from point sources, we have varied only the parameters of the 
source J0334. 3-3728 very close to the Fornax cluster as nuisance parameters. Pa- 
rameters of other point sources in the 2FGL catalog have been fixed to the best 
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Figure 2: Fitting results and limits from 7-ray observations. See text for explanation. 



fit values under background only hypothesis. With use of make2FGLxml .py [IB], we 
have included 2FGL sources within a radius of fifteen degrees as the background. 

We then put a limit through a signal plus background fit. For the experimental 
data, we have used P7S0URCE_V6 event class of the Fermi-LAT data measured be- 
tween August 4, 2008 and April 21, 2013 in the energy region E > 400 MeV, and 
in the region of interest (ROI) with a radius of ten degrees centered at the clus- 
ter. We have selected events passing the conditions recommended by the Fermi-LAT 
collaboration, using the package Fermi Science Tools [TT] . 

3.2.4 7-ray from the wino annihilation 

We have also considered the limit on the wino mass from the wino annihilation. 
This is because that the annihilation cross section of the wino (mainly into W boson 
pair) is significantly enhanced in some 0{1) TeV mass region. The most severe 
limit is again from several gamma-ray observations at Fermi-LAT. We have therefore 
included the limits in our analysis obtained from observations of dwarf spheroidal 
galaxies [HI HH] , galactic clusters [l9l HH] , and diffuse photons [361 EI] • 
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3.3 Results 



We are now at the position to discuss out results on the electron/positron flux and 
limits from gamma-ray observations, all of which have been obtained from numerical 
computations based on the strategy mentioned in the previous subsection. In figureE], 
the fitting result is depicted in each panel as contour lines of Ax^ = 2.3, 6.0, and 9.2, 
corresponding to 68th, 95th, and 99th percentile of the Chi-Squared distribution with 
two degrees of freedom, respectively. The red, green, and blue contours are the results 
of fitting to the AMS-02, Fermi-LAT, and both the AMS-02 & Fermi-LAT data. The 
PAMELA data of the electron flux is always used in the all fittings, where it plays 
an important role to fix the normalization of the background electron flux. The 
limits from gamma-ray observations at the Fermi-LAT experiment [3S1 ED] are also 
shown in each panel, where the black solid line named "Fermi 7 (Isotropic)" and 
the blue dash-dotted line named "Fermi 7 (A11-GP-|-GC)" are from the isotropic 
component and all sky survey of the diffuse gamma-ray observations, while the pink 
dotted line named "Fermi 7 (Fornax)" is from the gamma-ray observation of the 
Fornax cluster. The limit form the wino annihilation is shown as a grey shaded 
region named "Fermi 7 (Annihilation)". The best fitted mass and width of the 
wino are summarized in table [1] with the absolute value of the divided by the 
degree of freedom. For convenience to estimate the uncertainty associated with 
electron and positron propagations in our galaxy, the fitting results using Ml and 
M2 diffusion models are also shown in the table. One can see that the uncertainty 
of the propagation model does not change the result drastically. In figure [3l we have 
also plotted the positron fraction and the combined electron plus positron flux for 
the case of best fit parameters with keeping the wino mass at 3 TeV. 

For the results, we first notice that the decaying wino dark matter with LLE'^ 
interactions is indeed consistent with cosmic-ray electron and positron fluxes reported 
by AMS-02 and other collaborations. Furthermore, the most of the parameter space 
(the mass and width of the wino) explaining the positron anomaly is still consistent 
with gamma-ray observations, which is sharp contrast to the case of the decaying 
dark matter with LHu interactions. It is also clearly seen that the wino dark matter 
decaying into leptons of first and second generations requires the shorter lifetime 
than the one decaying into third generation leptons as expected. On the other hand, 
in some cases, the limits from the diffuse gamma-ray observations start excluding 
the parameter region. When the observation of its isotropic component is (officially) 
updated using the full data set that the Fermi-LAT collaboration already has, the 
decaying wino dark matter with some particular LLE'^ interactions such as L3L2-EI 
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Figure 3: Some examples of cosmic-ray signals with the wino mass of rriDM = 3 TeV: Left 
panels show the predicted positron fraction and experimental data {U O [U 152] /. while 
right panels show the electron and/or positron flux f33t [351 E21 EH I55f . The shaded region 
shows the systematics error associated with the scale of absolute energy at the H.E.S.S. 
experiment. The best Etted Fermi energy scale is taken in the plot. 
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Table 1: Fit results in some selected parameter space. See text for explanation. 



Mode 


Data 


m^ino [GeV] 


logio(r [sec]) 


xVdof 


Background 


AMS-02 






55.1/67 


Background 


Fermi 


_ 


_ 


47.8/43 


Background 


AMS-02 + Fermi 


- 


- 


271.8/93 




AMS-02 


806 


26.8 


43.7/65 




Fermi 


2611 


26.3 


26.3/41 




AMS-02 + Fermi 


2290 


26.4 


65.9/91 




AMS-02 + Fermi 


3000 (fixed) 


26.4 


69.2/92 


L3L2E2 


AMS-02 


1581 


26.5 


40.3/65 


L3L2E2 


Fermi 


5055 


25.9 


12.3/41 


L3L2E2 


AMS-02 + Fermi 


4209 


26.1 


62.4/91 


L3L2E2 


AMS-02 + Fermi 


3000 (fixed) 


26.1 


70.3/92 




AMS-02 


2784 


26.0 


36.1/65 


-^3-^2-E'3 


Fermi 


7439 


25.5 


25.5/41 


L3L2-E3 


AMS-02 + Fermi 


7957 


25.8 


57.2/91 


L3L2-E3 


AMS-02 + Fermi 


3000 (fixed) 


25.8 


89.7/92 


L3L2^f(Ml) 

^ 1 \ / 


AMS-02 + Fermi 


2276 


26.3 


65.7/91 


L3L2^2'(M1) 


AMS-02 + Fermi 


4172 


26.0 


61.6/91 


L3L2^|(M1) 


AMS-02 + Fermi 


7887 


25.8 


56.8/91 


L3L2Ef{M2) 


AMS-02 + Fermi 


2322 


26.4 


74.5/91 


L3L2^|(M2) 


AMS-02 + Fermi 


4979 


26.1 


72.1/91 


L3L2^|(M2) 


AMS-02 + Fermi 


5678 


25.8 


63.8/91 



will be favored or excluded. Even if the update result does not show the deviation 
from astrophysical backgrounds, the decaying wino dark matter into first and/or 
second generation leptons, such as the decay via L1L2E2, will survive. 

Another interesting thing we should point out is that the mass of the wino dark 
matter can be estimated from the results. It is trivial that the wino mass much 
lighter than 1 TeV is not consistent with the observed positron fraction. On the 
other hand, the wino mass much heavier than 1 TeV is also disfavored because the 
lifetime is too short to evade the gamma-ray constraints. Roughly speaking, the wino 
mass should be at around 3 TeV, which is nothing but the region the pure gravity 
mediation model predicts. Furthermore, with an appropriate LLE'^ interaction, the 
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best fit value of tlie wino mass is just around 2.7 TeV, wliicli is a very interesting 
coincidence because tlie pure wino dark matter predicts tlie mass of 2.7 TeV wlien 
its abundance observed today is from its thermal relic. 

Finally, we would like to comment on a slight discrepancy between the result 
fitting to AMS-02 and that fitting to Fermi-LAT. The discrepancy is actually not 
due to our model but from the one between the two experimental results as already 
pointed out by recent studies [56]. The essential reason is that the positron fraction 
that the AMS-02 collaboration reported is somewhat shallower than that reported 
by the PAMELA collaboration which is known to be consistent with the combined 
electron plus positron flux reported by the Fermi-LAT collaboration. As a result, 
for the decaying wino dark matter scenario, the wino lifetime to explain the AMS-02 
result is always longer than that to explain the Fermi-LAT result, and it makes the 
discrepancy. The AMS-02 experiment will report on their new results of electron 
and positron fluxes in near future, and then we will have a clue to understand the 
discrepancy. Once we know the reason, more accurate analyses of the decaying wino 
dark matter will be possible. 

4 Summary 

In light of the recent result of the positron fraction reported by the AMS-02 col- 
laboration, we have discussed the decay of the wino dark matter through LLE'^ 
interactions in the framework of the PGM model. The PGM model predicts the 
almost pure neutral wino as a lightest supersymmetric particle, namely a candidate 
for dark matter, whose mass is predicted to be around 1 TeV. This result is very 
interesting because the AMS-02 result in fact requires the dark matter mass around 
1 TeV if the excess of the positron fraction is interpreted as a decaying dark matter 
signal. Furthermore, we have explicitly shown that the decaying wino through LLE'^ 
interactions well reproduce the AMS-02 result with being consistent with constraints 
from anti-proton and gamma-ray observations. We have also discussed a model of 
the R-parity violation where only leptonic R-parity violating interactions are gener- 
ated in the appendix, while neither the other hadronic R-parity violating interactions 
nor bi-linear R-parity violating interactions are generated. Such a model is crucial 
to avoid the constraint from cosmic-ray anti-proton observations. 

Here, we comment on some future prospects for testing the decaying wino dark 
matter. One of the distinct properties of the wino dark matter is that both its 
annihilation and decay provide signals to cosmic-ray and gamma-ray observations. 
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This is sharp contrast to other dark matter candidates such as the decaying gravitino 
or stable particles whose decay are forbidden. The decaying wino dark matter can 
therefore be tested in various indirect searches. For instance, the wino produces anti- 
protons thorough its annihilation, which is expected to be detected at the AMS-02 
experiment in near future O [9] . As can be seen in previous section, the wino also 
produces diffused gamma-rays mainly through its decay, which can also be used to 
test the dark matter when data is accumulated enough. 

Another interesting test for the decaying wino dark matter is the observations of 
galactic clusters and dwarf spheroidal galaxies using future gamma-ray telescopes. 
The decay process of the wino dominantly contributes to the signal from galactic 
clusters, which will be detected in the future at the next generation air-Cherenkov 
telescopes [21]. On the other hand, the annihilation process of the wino becomes 
comparable to the decay process at the signal from dwarf spheroidal galaxies, which 
is also expected to be observed in future satellite experiments [5?]. When both 
signals are detected and their energy spectra are measured accurately, these signals 
from the above two processes have a potential to be a smoking-gun signature. 
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A Possible origin of LLE^ interactions 

In this paper, we have assumed that the R-parity violation is dominated by LLE'^ 
interactions. If we consider the GUT, however, LLE'^ interactions are accompanied 
by the other R-parity violating interactions U'^D^D'^ and QD'^L when they originate 
in the R-parity violating interaction 10 5* 5*. Here, we have used SU(5) GUT repre- 
sentations for the MSSM matter fields, 10 = (Ql, Er) and 5* = {Dr, Ll). Note 
that those accompanied R-parity violating interactions do not cause rapid proton de- 
cay because of A = 0(10^^^) which is favored for the observed excess of the positron 
However, they cause a large flux of cosmic-ray anti-protons which easily 



fraction 
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"^^This should be contrasted to the decaying gravitino dark matter scenario [S51[3T] which explains 
the excess of the positron fraction through the R-parity violating interactions. That is, since the 
decaying gravitino dark matter requires much larger R-parity violating LLE'^ interactions, the 
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10 


5* 


$(5) 


${5*) 




X{1) 


S(l) 










SU(2)j^ 


1 


1 


2 


2 


2 


2 


3 + 1 


1 


1 


2 


2 


U(1)h 








-1/2 


1/2 


-1/2 


1/2 





1 


-1 


-1/2 


1/2 



Table 2: Charge assignments on the SU(5) x U(2)h GUT gauge groups. Flavor indices 
of 10, 5*, Lh, Lh, Eh, and Eh are suppressed for simplicity. 

exceeds the observed constraint on it [22]. 

In this appendix, we discuss a model where LLE'^ interactions are accompanied 
by neither U'^D'^D'^ nor QD^L interactions in a GUT-consistent way. For this pur- 
pose, we consider an SU(5) x 11(2)// product group GUT model [SH], where LLE'^ 
interactions originate not in 10 5* 5* but in R-parity violating interactions consisting 
of additional lepton-like fields, Lh, Lh, Eh and Eh, in the \J{2)h sector (see table 
2)0 

Let us first briefiy review the product group GUT model based on SU(5) x \J(2)h- 
In this model, no adjoint of SU(5) is introduced and the GUT gauge symmetries are 
broken by the vacuum expectation values of bi-fundamental fields $ and $ in table [2], 
($^) = v6l and (<l'^) = v6^, where v denotes a dimensionful parameter around the 
GUT scale and indices run as a = 1-5 and i = 1-2. Such vacuum expectation values 
are obtained as a supersymmetric solution of the superpotential: 

W = V2^] (t^)| $1 + V2Xj (t^)| x'-V2v^j:°, (4) 

where denotes Pauli matrices divided by two (t^'^'^ = cr^'^'^/2) and the two- 
dimensional unit matrix divided by 2 (t° = I2X2/2). With the vacuum expectation 
values, SU(5) x \J{2)h are spontaneously broken down into MSSM gauge groups. 

As a remarkable feature of this model, MSSM Higgs doublets are not embedded 
in the fundamental representations of SU(5) but embedded in the doublets of 1]{2)h, 
i.e. Hd ~ "I^X* and Hu ~ ^a^i- Therefore, the doublet-triplet problem is absent 
since no colored Higgs multiplets are introduced in this modelEj MSSM Yukawa 
interactions are effectively given by higher dimensional operators, 

W = — 105*$X + — 1010$X, (5) 

where fiavor and gauge indices are suppressed. These higher dimensional opera- 
tors may be generated by integrating out heavy particles, so that the dimensionful 
parameter M* is of the order of the GUT scale. 



accompanied U^D^D'^ and QD'^L interactions lead to a too short lifetime of a proton. 
*^For a related model of R-parity violation based on a product group GUT, see reference 
'^^AU the other modes consisting of $'s and S obtain masses at the GUT scale. 
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Let us now introduce R-parity violation consisting of the hidden lepton-hke fields, 



W = XhLhLhEh, (6) 

where couplings Xh is a dimensionless parameter. The hidden lepton-like fields have 
masses Ml and Me which are assumed at GUT scale. Note that we need at least two 
pairs of the lepton-like doublets for non-vanishing R-parity violating interactions. 

After the GUT breaking, the R-parity violation of the lepton-like fields is medi- 
ated to the MSSM leptons through W = Lk Lh, which cause mixing between 
Lh and L's with the angles of Hkv/Mi. In addition, we also consider higher dimen- 
sional operators, W = {k'^/M^) 10^ = {k'i^v'^ /M^)E^ Eh, which lead to the 
mixing angles between Eh and £"s, k^v"^ /{M^Me)- Through these mixings, R-parity 
violating interactions of the lepton-like fields generate the required leptonic R-parity 
violating operator, 

Xh 4 
^ - MIMeM^ ^'^^^'^ 

while U^D^D'^ and QD'^L interactions are suppressed. In this way, we obtain a model 
in which only the leptonic R-parity violation is realized in a GUT consistent manner. 
For Kfc = Ml ^ Me V ^ 10^^ GeV and 10^^ GeV, and A^^ ~ 10-^^ we 

obtain the required R-parity violation A ^ 10~^^, 
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